Fermionic Stochastic Schrodinger Equation and Master Equation: An Open System 
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This paper considers the extension of the non-Markovian stochastic approach for quantum open 
systems strongly coupled to a fermionic bath, to the models in which the system operators commute 
with the fermion bath. This technique can also be a useful tool for studying open quantum systems 
coupled to a spin-chain environment, which can be further transformed into an effective fermionic 
bath. We derive an exact stochastic Schrodinger equation (SSE), called fermionic quantum state 
diffusion (QSD) equation, from the first principle by using the fermionic coherent state representa- 
tion. The reduced density operator for the open system can be recovered from the average of the 
solutions to the QSD equation over the Grassmann-type noise. By employing the exact fermionic 
QSD equation, we can derive the corresponding exact master equation. The power of our approach 
is illustrated by the applications of our stochastic approach to several models of interest including 
the one-qubit dissipative model, the coupled two-qubit dissipative model, the quantum Brownian 
motion model and the N-fermion model coupled to a fermionic bath. Different effects caused by the 
fermionic and bosonic baths on the dynamics of open systems are also discussed. 

PACS numbers: 03.65.Yz, 42.50.Lc, 05.40.-a 



I. INTRODUCTION 



The theory of open quantum systems has experienced 
a resurgent interest because of the rapid development 
of quantum experimental technologies and their applica- 
tions to the fabrication and manipulation of quantum 
devices (e.g. photonic devices, quantum dots, nano- 
mechanical oscillators). However, an intricate problem 
exists since in reality no system can be completely iso- 
lated from its environment (bath, reservoir etc), and the 
dynamics of the system of interest will be profoundly af- 
fected by the couplings to its environment 0, 0| . When 
the quantum open systems are coupled to a Markov envi- 
ronment, the Lindblad master equation is a critical tool 
which can be used to study the dynamics of the open 
systems [3(. When the Born-Markov approximation is 
no longer valid, namely, the coupling between system 
and environment is not weak and the environment can- 
not be approximated by a broadband bath, one must 
extend the standard Markov theory to a more general 
non-Markovian environment. Several attempts to derive 
the evolution equation of open quantum systems beyond 
the Markov approximation have been proposed @, 0- 
8] . Notably, the non-Markovian quantum state diffusion 
(QSD) approach developed by Strunz and his cowork- 
ers has showed momentous potential of solving large sys- 
tems (multi-qubit or multi-cavity) [9h15j]. Moreover, as a 
computing tool, many numerical advantages of the QSD 
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approach permit its use in several domains such as high- 
precision measurement [161 ]. entanglement dynamics [171 ] 
and coherence dynamics of the large molecules in bio- 
physics [18| etc. Therefore, it is highly desirable to extend 
the QSD approach for the bosonic baths to the fermionic 
case where the non-Markovian features have played in- 
creasingly important role [19l - 423j | . 

The primary theme of our current paper is to establish 
an exact quantum approach for a class of quantum sys- 
tems interacting with a fermionic bath. We will consider 
a class of systems such that the systems and fermionic 
bath are distinguishable, hence the system Hamiltonian 
and the bath operators commute. The system of interest 
in this case may consist of one or more effective particles 
such as spins, effective fcrmions etc. The case where the 
system and bath operators ant-commute will be covered 
in-depth in a separate paper [2J]. It is noted that the 
commutative model we proposed arises from many phys- 
ical settings including spin bath and fermionic bath (e.g. 
see Appendix A). 

We will derive a fermionic stochastic Schrodinger 
equation for an open quantum system embedded in a 
fermionic bath, called the fermionic QSD equation. To 
illustrate the power of our approach, we solve several 
models as examples by using this new technique, includ- 
ing a one-qubit dissipative model, a two-qubit dissipative 
model, the quantum Brownian motion in a fermionic bath 
and a multiple-particle model. In the first example, we 
give the explicit analytical solution without any approx- 
imation in a special case. In the second example, we 
show how to construct the crucial Q operator contained 
in the fermionic QSD equation. In the third example, 
we consider a continuous variable model where a Brow- 
nian particle is immersed in a bath of fermion particles. 



The last example involves a genuine multi-particle sys- 
tem that has been solved exactly by our QSD approach. 
Finally, the difference between the bosonic bath and the 
fermonic bath is discussed. 

This paper is organized as follows. In section [Til we 
introduce the general commutative fermionic bath model 
and derive the fundamental dynamic equation for this 
model. In section Mil we derive the formal exact mas- 
ter equation from the QSD equation. In section IIV1 we 
present a simple example of using this fermionic QSD 
approach to solve the one-qubit dissipative model. In 
section [V] we solve the two-qubit dissipative model to 
show the construction of some complicated Q operator. 
In section IVI| we apply our fermionic QSD approach to 
a continuous variable model to solve the quantum Brow- 
nian motion model in a fermionic bath. In section I VIII 
we solve a genuine multi-partite system, the N-fermion 
model, to show that our approach is not only applicable 
to small systems (one-qubit or two-qubit), but is also ap- 
plicable to large quantum open systems. Based on the 
last example, we also evaluate the differences between 
the bosonic bath and the fermonic bath. Finally, in sec- 
tion lVHll we conclude the paper. In Appendix A, we pro- 
vide an effective commutative model consisting of spin- 
less fermions as an environment. In Appendix B and C, 
we present the details of derivation of the non-Markovian 
QSD equation for a fermionic bath. In Appendix D, we 
prove a Novikov-type theorem for a Grassmann Gaus- 
sian stochastic process, which plays a crucial role in de- 
riving the exact master equation from the corresponding 
stochastic Schrodinger equation. In Appendix E-F, we 
derive explicit equations of motion for the coefficients of 
the master equation for the examples presented in this 
paper. 



II. NON-MARKOVIAN QSD EQUATION FOR 

AN OPEN SYSTEM COUPLED TO A 

FERMIONIC BATH 

For a quantum open system interacting with a 
fermionic environment, the total Hamiltonian may be 
written as 



fft, 



H a + Hy 



H„ 



(1) 



where H s is the Hamiltonian of the system, Hy, is the 
Hamiltonian of the bath and H- mt is the interaction term. 
When we consider the fermionic bath, Hy, and Hint can 
be written as (setting h — 1 throughout the paper) 



H h = y^Ujc\ci, 



H int = Y^(aZ4L + 9iHci), 



(2) 



(3) 



where c\ and Ci are fermionic creation and annihilation 
operators {£i, cA = Sij. Here we emphasize that the bath 



may consist of a set fermions or spins (e.g. see [26|; An 
example is shown in Appendix A.) 

In the interaction picture, the total Hamiltonian be- 
comes 

H tot (t) = H S + ^2(g*e iUit c]L + 9i e- iUit L%). (4) 



We use the fermionic coherent state (e.g., see [27|, [2g|) to 
describe the state of environment. For a single mode, the 
fermionic coherent state is defined as 



q|&) = &!&>, 



(5) 



where £j is a Grassmann variable which satisfies the fol- 
lowing properties {&, £,-} = 0, {£i, £*} = 0. Generally, the 
coherent state can be expanded in terms of Fock states 
as l£») = |0) — Ci c l |0). The coherent states (ket and bra 
vectors) for the multi-mode environment are given by 

I0 = |6)®|6)®l&>®...and<el = (ei|®(6|®(&l®.... 
Now, we can define 



MO = »«*(*)>, 



(6) 



where \iptot(t)) is the total state vector for the system and 
environment, and (£| is a coherent state representation 
for the environment. In this paper, we focus on the case 
where the initial state of the bath is vacuum state. The 
finite temperature bath will be discussed in [2J| by using 
the Bogoliubov transformation [29(. With the coherent 
state representation, we can derive the non-Markovian 
QSD equation for the fermionic bath as 



d_ 

at 



Mn 



[~iH s + LQ -til dsK(t, 
Jo 






where Q 



iJ2i9ie iUit ^ and K(t,s) 



(7) 



5TJJ<?i| 2 e lLJi ^ s ) is the correlation function. (De- 
tails of the derivation can be found in Appendix B) Wc 
use jpript{£,*) to denote the left-functional-derivative 
with respect to the Grassmann variables. Our fermionic 
approach is applicable to arbitrary correlation functions 
especially for the general non-Markovian case. 

Similar to the formal bosonic QSD equation [9], the 
fermionic QSD contains a time-nonlocal Grassmann func- 
tional derivative which renders a direct application of the 
derived fermionic QSD extremely difficult, if not impossi- 
ble. In order to find a time-local QSD equation, one can 
introduce a time-dependent operator (also ^-dependent 
in general) Q, defined as 






Q(t,s,nM? 



(8) 



If no confusion arise, we will use the shorthand notation: 
Q = Q(t,s,£*). With this Q operator, the exact stochas- 
tic QSD equation can be written as 



_5 

at 



Mn 



iH s + l^-Dq]MC), 



(9) 



where Q{t,£*) = J Q dsK(t,s)Q(t,s,£*). The stochastic 
QSD equation for a fermionic bath we have presented 
here is an exact equation of motion for the open quan- 
tum system directly derived from the microscopic Hamil- 
tonian without any approximation. It should be noted 
that in our derivation of the QSD equation, we have not 
explicitly specified the system Hamiltonian and the cou- 
pling operators, H s and L. Here we have introduced a 
new type of stochastic process £ t *. The solution of our 
QSD equation is called a Grassmann quantum trajectory. 
By construction, the reduced density matrix of the open 
system can be recovered by the statistical mean over the 
Grassmann noise. Although the fermionic QSD equation 
looks formally similar to the bosonic case, the dynamic 
behaviors of the system governed by the two types of 
equations can be different due to distinct differences be- 
tween the bosonic and fermionic particles. Mathemat- 
ically, the most striking difference between the bosonic 
and fermionic QSD equations is that the former contains 
a complex Gaussian noise while the latter is driven by 
a non-commutative Grassmann Gausian noise. We will il- 
lustrate the difference in a concrete example in a subse- 
quent section. 

In order to derive the dynamic equation for the Q op- 
erator, we consider the consistency condition (CC), 



fact, we may expand Q operator as 



A! 

<*£ dt 



Mn 



1A 

dt 8$ 



mn- 



(10) 



Applying the QSD Eq. (|9|) to CC, the equation for Q 
operator is derived as 



di 



Q = -i[H B ,Q]-{LQ,Q} 

-frQ{-e)Q + QpQ-fr-AQ, (11) 



where the sign of Q{— £*) depends on the functional form 
of noise contained in Q. (Details of derivation and discus- 
sion can be found in Appendix C.) The initial condition 
for the Q operator is 



Q(t,s = t,t)=L. 



(12) 



However, for a simple case that Q is independent of 
Grassmann noise, the equation for the Q reduces to 



^ t Q = -i[H s ,Q]-{khQ}-pQ,Q]. 



(13) 



Eq. (fTTj) and Eq. ((T3|) can be used to determine the ex- 
act Q operator. However, for most practical problems, it 
may be a daunting task to determine the exact Q. There- 
fore, it is important to develop a perturbation approach 
similar to that developed for the bosonic bath [lOf. In 



Q(M,n = Q (U, (M)+ / Q ti} (M,si)C ds i 



)Wl 
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Q (2) (t.*,*i,«2)£ 1 £ a daida 2 - 

o Jo 
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(n), 



JO 



'(*, s, si,...s n )£ 1 ...€* n d8i...ds n 

(14) 



Substituting this equation into Eq. (fTTj). one can derive 
the dynamic equations of the coefficients for each order 
Q( 1 \ Particularly, the zeroth-order term Q^(t, s) will 
satisfy the following equation (neglect all the noise terms) 

^Q i0) (t,s) = -i[H s ,QW(t,s)] - [Ltg(°)(t),Q(°)(i, s )], 

(15) 
where Q (0) {t) = f*Q ( -^(t,s)K(t,s)ds, and the initial 
condition is 



Q {0) {t,s = t) = L. 



(16) 



III. FORMAL EXACT MASTER EQUATION 

FOR AN OPEN QUANTUM SYSTEM COUPLED 

TO A FERMIONIC BATH 

Now, we will derive the master equation governing the 
reduced density operator of the open quantum system 
from the stochastic QSD equation ([9])- First, we define 
the stochastic density operator as 



p t = mn) (M-s)\ 



(17) 



It is easy to verify that the reduced density matrix of the 
open system can be written as 

P = ^2(n\Aot)(Aot\n) 

]J.dti*dS i e-^V^(n\0(Z\i>tot)('iptot\n) 

n 
J l n 



= (Pt. 



(18) 



where (...) s denotes the statistical mean over the Grass- 
mann Gaussian noise defined by 

(->. = yn^d&e-^H..). (19) 

From this expression, we say that the reduced density 
matrix can be unraveled by a set of Grassmann quantum 
trajectories |^t(£*))- 



From the QSD equation Eq. ©, we have 

d_ 
di 

thus, 



(M-0\ = (M-t)\ [iH". - ttft - QHt, -OL], (20) 



d -- d lf>\ 

= ((-iH s + Le t -^Q)P t ) s 

= -i[H s , p] + L(QP t ) s - (P t £t).tt 
-V{QP t ).-{P t QH-Z)).t. (21) 

In order to establish the exact master equation from 
the fermionic QSD equation ©, one needs to handle 
the terms (Pt£,t)s etc. In fact, we can prove a Novikov- 
type theorem for the Grassmann Gaussian noise (see Ap- 
pendix D), 



(P&)s = -(QPt)., 

(GP t ), = (PtQH-O), 



(22) 
(23) 



With the help of the Novikov-type theorem for the Grass- 
mann noise, the exact master equation can be written as 

J^p = -i[H„ p] + [L, (PtQHS)).] + [{QPt).,ft]. (24) 

If the operator Q is independent of the Grassmann noise, 
then the exact master equation is immediately obtained, 



d_ 
di 



p = -i[H s ,p] + [L,pQi} + [Qp,fr] 



(25) 



Moreover, in the Markov limit, Q — jfL, this master 
equation reduces to the standard Lindblad master equa- 
tion: 

^ t P=-i[H s ,p}+ 7f [L,pfr}+ 7f [Lp,P}. (26) 

In subsequent sections, we will derive several interesting 
master equations from the corresponding QSD equations. 



IV. EXAMPLE 1: ONE-QUBIT DISSIPATIVE 
MODEL 



H s = -a z 



H h = y^Ujclcj 



H int =Y,^tc\L + g i Dc i ), 



4 

(28) 

(29) 
(30) 



where L — &- for this particular model. 

From Eq. ([13]), the solution for the Q can be obtained 
as 

Q(t,s)=x 1 (t,s)a-, (31) 

with the initial condition 

Q{t,s = t) = L = fr-, (32) 

and the coefficient x\{t, s) is shown to satisfy 

— xi(t, s) = [iu + Xi(t)]xi(t, s), (33) 

at 

where Xi(t) — Lxi(t, s)K(t, s)ds, and K(t, s) is the cor- 
relation function, and the initial condition is given by 
X\(t, s) = 1. 

Thus, the exact Q operator can be fully determined. It 
is worth noting that this Q operator has the same form as 
the bosonic case [lOj. Finally, the explicit QSD equation 
for this model is 

§- t MC) = [-*|*, + &-% - Xi(«)a + a_]V*(0> ( 34 ) 
and the exact master equation is 
d 



dr 



= -i[H s ,p} + [L,pQi} + [Qp,fr] 

= -i-ziZzP - P^z) + X*(t)(a-p& + - pcr+cr-) 



X 1 {i)(a-pa + - a+a-p). 



(35) 



With this exact master equation, the dynamics of this 
model can be fully determined. 



We start with a very simple example, one-qubit in 
fermionic bath. This is a special case where it is pos- 
sible to derive the fully analytical solution without any 
approximation. 



A. Master equation and non-Markovian quantum 
dynamics 

The total Hamiltonian for the one-qubit dissipative 
model may be written as 



B. A limiting case — the environment consists of 
only one fermion 

Now, we consider a very special case for the one-qubit 
model where the "environment" [3Q] contains only one 
fermion. By analytically solving this model, we show 
explicitly that the fermionic QSD gives rise to identical 
results to those predicted by the ordinary quantum me- 
chanics. The model is described by the following Hamil- 
tonian, 



Htnt — H s + Hh + H„ 



(27) 



H, 



-a z + ui b c f c + (g*(T_c t + gcr+c), (36) 



and the zero-temperature correlation function becomes 



K(t,s) = \g\ 2 e' iu ">^- S \ 



(37) 



Substituting the correlation function into the expression 
of Xi(t) — Lxi(t,s)K(t,s)ds, we will find the differen- 
tial equation for X\ (t) as 

9 X 1 (t) = \g\ 2 - WiXi.it) + iuX x {t) + X^t) 2 . (38) 



Of 

For simplicity, we consider the resonance case, then the 
solution X\ (t) can reduce to 



X 1 (t) = \g\tan(\g\t). 



(39) 



From the master equation Eq. Q35p . we can calculate 
time evolution for the off-diagonal elements in the density 
matrix. 

|p 2 i = | <*+) = Tr{^pa + ) = iup 21 Xl(t) hl . 

Finally, we can derive the time evolution for p 2 i as 

M) = fa(p)e**co8[\g\t]. (41) 

Similarly, we can get 

Pi2(t)=/Sia(0)e- fart cos[| fl f|t]. (42) 

This result shows that the coherence (off-diagonal ele- 
ments in density matrix) will decrease and increase peri- 
odically. 

On the other hand, we can easily solve this simple 
case using elementary quantum mechanics. Since this 
is only a two-body problem, we can solve the evolution 
for the whole system in a straightforward manner. One 
can check that elementary quantum mechanics gives rise 
to the identical results obtained by the fermionic QSD 
approach in Eq. (J41I42I) . 



two-level models) interacting with a common fermionic 
bath. We will show how to construct exact and approxi- 
mate Q operator in this example. The total Hamiltonian 
of this model can be written as, 



Htot — H s + Hb + H-, 



int. j 



(43) 



where 



h s = oj a &^ + Lose? + J xy {o+v B + a A af) + J z v?fii 



flint = £(fljc}£ + gjCjtt), 



(44) 



Here, the operator L = kao^ + kb&- describles the pat- 
tern of interaction to the environment, ka and kb are 
constants describing different coupling strengths for the 
two qubits. 

The perturbative zeroth-order Q operator can be de- 
rived as 



Q (0) (M) = ^/,(t, s )& 



(45) 



i=l 



where Qi (i — 1,2,3,4) are the time-independent basis 
operators, and fi(t,s) are time-dependent coefficients. 

The four basis operators in terms of the Pauli matrices 
may be written as 



V. EXAMPLE 2: COUPLED TWO-QUBIT 
DISSIPATIVE MODEL 



Oi = of, Q-2 = o u , Q 3 = df&Z, Qa = &?6±, (46) 



In this section, we consider a system containing a pair 
of coupled two- level systems (spins or some other effective 



J 



From Eq. (|15p . we can derive the differential equation 
for the coefficients as 



|/ 2 (M) = 
|/a(M) = 



+2iu A fi - Uxyh + 1iJ z fa + K A Fifi - kb-Fi/s + K B F 3 fi + n B F 3 f 4 + K B F 4 f 3 + k a FaU, (47) 

+2iuj B f2 - iJxyh + 2iJ z f 3 + n B F 2 f 2 - n A F 2 f 4 + n B F 3 f 3 + n A F 3 f 4 + n A F 4 f 2 + n A F 4 f 3 , (48) 

+2iu B f 3 - iJ xy h + 2iJ z f 2 - K A F 2 fi + n B F 2 f 3 + K A F 3 h + K A F 4 f 2 + K B F 3 f 2 + n A F 4 f 3 , (49) 

+2iu A f 4 - iJxyf-2 + 2iJ z fi - n B Fif 2 + K A Fif 4 + n B F 3 fi + n B F 3 f 4 + n A F 4 f 1 + n B F 4 f 2 , (50) 



where F^t) = f* dsK(t,s)fi(t,s) (i = 1,2,3,4), and the 
initial conditions are 



fl(t,S = t) = KA, 
fs(t,S = t) = K B , 
f 3 (t,S = t) = 0, 

f 4 (t,s = t) = 0, 



(51) 
(52) 
(53) 
(54) 



Moreover, we can also determine the exact Q operator 
for this two-qubit model. We can verify rigorously that 
the exact Q operator contains five terms where the last 
term is noise-dependent. The details of the derivation is 
presented in Appendix E. If we use this zeroth-order Q 
operator, the master equation can be explicitly written 
in the following form 



dt 



p = -i[H s p- pH s ] 

4 

+ iY, f:\lpq\- 



pom+H.c.}. 



(55) 



Next, we consider a simple case, in which all the param- 
eters are symmetric for two qubits, i.e. lua = u>b = w, 
ha = hb = 1. Then, we can derive the following master 
equation 



dt' 



-i[H Bt p] + [L t pQi] + [Qp,ft] 



sf)p- 



■p{K 



a B a A )p- 






piaiai 



&»<jA)] 



U xy [(cr+cr B 

iJ z [<J z <j z p — pa'^a 

{F*[(<y A + a B )pai-pai(a A +a B )} 

F*[(a A + a B )pa B -pa B (* A + a B )] 

F*[(a A + o B )po A ol pa A * B (a A + a B )} 



m&t 



o B )po z o+ 



A --B~As~A 

pa z a + (a_ 



a B )}+H.C.}. 
(56) 



The master equation derived above is valid for a gen- 
eral correlation function. For numerical simulations, one 
need to consider a special example of the correlation func- 
tion. A general correlation function may be written as 



K(t, s) = / du J(uj)[coth(u/2k B T) cosu(t - s) 
Jo 



— i sinw(i — s)] 



(57) 



where J(u>) is the spectral density. If we choose J(u>) — 



Tue 



-(■?-) 



which is so-called Ohmic case, the correlation 



function in the zero-temperature can be written as 

r 



K(t,s) 



+ i(t-s)}< 



(58) 



where uj c is the cut-off frequency. 
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FIG. 1: Time evolution of concurrence for different ui c . 
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FIG. 2: Time evolution of |Fi(£)|, \F 2 {t)\, \F 3 (t)\, and |F 4 (i)|. 
In the symmetric case, \Fi(t)\ - \F 2 (t)\, \F 3 (t)\ = |F 4 (i)|. 
The parameters are Wi = CJ2 = W = 1, 



VI. EXAMPLE 3: QUANTUM BROWNIAN 
PARTICLE IN A FERMIONIC BATH 



We consider a continuous model consisting of a Brow- 
nian particle interacting with a fermionic bath. The 
Hamiltonian of the Brownian particle is given by, 

H s = u m {p 2 + q 2 ). (59) 

The Hamiltonian of the fermionic bath is 

H h = ^Wifijci, (60) 

i 

and the interaction Hamiltonian is given by 

Hint = q^2{g*cj + 9iCi)' (61) 



So, the total Hamiltonian is 

H to t = H s + Hb + #int ■ 



(62) 



Applying our QSD approach to this model, it can be 
easily shown that Q operator takes the following form: 

Q = xi(t, s)q + x 2 (t, s)p + x 3 (t, s,C)pq 

+ Xi{t, s, C)P 2 + x B (t, s, Of + ■■■• (63) 

which is a infinite scries, therefore, it is difficult to de- 
termine the exact Q operator. A useful approximation is 
to neglect all the noise-dependent terms, after which we 
obtain the so-called zeroth-order approximate Q as 



Q « xi(t, s)q + x 2 (t, s)p. 



(64) 



Substituting this approximate Q operator into Eq. (|15|) . 
we can derive the differential equations for the coefficients 
X\(t, s) and x 2 (t, s) as 



d_ 
at 

at 



xi(t,s) = 2uj m x 2 (t,s) +iX 2 (t)xi(t,s) - 2iXi{t)x 2 {t,s), 

(65) 



x 2 (t, s) = —2u m Xi(t, s) - iX 2 (t)x 2 (t, s). 



(66) 



The initial conditions for coefficients Xi(t, s) and x 2 {t, s) 
are 



X\(t, s = t) = I , 
x 2 (t,s = t) =0. 



(67) 
(68) 



Using this approximate Q operator, the master equation 
can be written as 

j t P = -i[H a ,p] + [L,pQt} + [Qp,fr] 

= -ioj m [(p 2 + q 2 )p- p(p 2 + q 2 )} 

+ {X{ [qpq - pqq] + X* [qpp - ppq] + H.C.}. (69) 

It should be noted that the A| (including its conjuga- 
tion X 2 ) does not exist in the Markov limit. Hence, 
the approximate Q operator defined above is different 
from the Markov approximation. It is also different from 
the weak-coupling approximation since the approximate 
Q still contains the higher-order terms of the coupling 
constant. We expect that the master equation obtained 
from the approximate Q will be valid in a weakly non- 
Markovian regime. 

From the master equation we can derive the evolution 
equations for all the mean values of operators q 7 p, 



— (q) =2w m (p), 

d 

— {p) = -2u m (q) 



iXt{q)-iX* 2 ($)+iX x {q) 



(70) 

iX 2 {p). 
(71) 




FIG. 3: Time evolution of mean values of operator q in dif- 
ferent environments. The parameter 7 indicates the memory 
effect. The other parameters are uj m — OJ = 1, Q = n/2. 



In Fig. [3j we plot the time evolution of (q) in dif- 
ferent kinds of environments with different 7. In order 
to show the transition from non-Markovian to Marko- 
vian, Ornstein-Uhlenbeck noise K(t,s) = 2 e ~( 7+i ^l*~ s l 
is chosen in our numerical simulations. The reason of 
using Ornstein-Uhlenbeck noise is that the memory time 
of the environment can be described by one parameter 
I/7. Fig. [3] shows how the evolution of (q) is affected by 
7. This is a unique phenomenon in the non-Markovian 
case. 



VII. EXAMPLE 4. N-FERMION SYSTEM 
COUPLED TO A FERMIONIC BATH 



Dynamic equation for the general N-fermion 
model 



In the last example, we will establish the exact time- 
local fermionic QSD equation and master equation for a 
genuine multipartite system coupled to a fermionic bath. 
We show that, using the fermionic QSD approach, the 
exact Q operator of the N-qubit model can be easily de- 
termined. 

More specifically, let us consider the following Hamil- 
tonian 



^tot — H g + H b + -ffint , 


(72) 


N b 




j—i 


(73) 


N b 





H int = Y,9 3 {c]L + Dc/). 
3=1 



(74) 



where H^ is a fermionic bath. We assume that the system 



of interest consists of N fermions, i.e. 

N, 



H s = y^ Ajalai 



(75) 



i=\ 



here, d\ and a% are also fermion creation and annihilation 
operators; the Lindblad operator is 

A 1 . 



^ = E 



a-. 



(76) 



This Hamiltonian could be an effective Hamiltonian 
transformed from a set of spins. For example, sup- 
pose that we have a long chain with N sites, if the first 
N S (N S < N) sites are treated as system and the other 
Nb sites are treated as bath (N s + Nb = N), then per- 
forming the Jordan- Wigner transformations for both the 
system and the bath, we may result in this type of effec- 
tive Hamiltonian (for details, see Appendix A). 

We can show that the exact Q operator of this model 
takes the following form 



N a 



Q = y^^Xi(t,s)ai, 



(77) 



and the differential equation for the coefficients in Q op- 
erator are given by 



0_ 
di 



N s 



Xj(t,s) = iAjXjjt, s) + y^ Xj(t)xj(t, s), (78) 



where Xj(t) = f K(t,s)xj(t,s)ds. So, Q(t) — 
J2i=i-X-i(t)ai. The exact master equation of this model 
is 

^- t P = -i[H s ,p} + [LJQi] + [Qp,P] 



Ns 



N a 



= -i[H 8 ,p] + [J2ai,P^X*(t)i 



N a 



A's 



[(x;^(*)6i)p.E a ^ 



(79) 



B. Fermionic versus bosonic baths 

It is instructive to consider a simple case with two 
fermions in the system (N s = 2). The Hamiltonian is 
then given by 

H s = bJia\ai + uj 2 a\d 2l (80) 

L = ai+a 2 , (81) 

it is easy to show that the exact Q operator is 

Q = X 1 {t)ai+X 2 {t)aa (82) 



where X\ (t, s) and X 2 (t, s) can be determined in Eq. (|78|) 
as TV;, = 2 case. Then, the explicit master equation can 
be written as 



dt 



p = — iu\{a\aip — pa\ai) — i^ia^^p — P&l^) 

+ {X*(t)(aipa\ - pa\ai) + Xl(t){a 2 pa\ - pa\a 2 ) 
+ X 2 (t)(aipa\ — pa\ai) + X 2 (t)(a 2 pa 2 — pa\a 2 ) 
+ H.C.}. (83) 



On the other hand, we can also solve this model ex- 
actly if the two effective fermions (spins) are coupled to 
a bosonic bath. The Hamiltonian takes the same form as 
Eq. (|72H76p . except that dj (ct) represent bosonic anni- 
hilation (creation) operators (also consider N s = 2 case). 
Using the non-Markovian QSD approach for bosonic bath 
|9| , the bosonic QSD equation can be derived as 

^tfcCO = [-*#. + lz *t - How*), (84) 
at 

where 0(t,z*) = J Q dsK (t, s)0(t, s, z*). In the bosonic 
QSD equation, the noise z\ = —i^ngXe^^z* is the com- 
plex (not Grassmann) Gaussian noise. The exact O op- 
erator is determined as follows 



0{t,z*)=X 1 {t)di+X 2 (t)a 2 +X ? ,{t)a\dia 2 

+ X4(t)a,2did 2 + i / ds'X^(t, s')z*, did 2 . (85) 



Details about the coefficients can be found in Appendix 
F. 

We use this particular example to illustrate different 
aspects between the fermionic and bosonic baths. As 
shown above, we can find the exact Q (O) operators for 
both the fermnionic bath and the bosonic counterpart. 
Since the exact dynamic evolution of the system will be 
fully determined by Q (O) operators, so we may com- 
pare the difference between the two operators given in 
Eq. ([82)1 and Eq. (|85|) . respectively. The first two terms 
Xi(t, s) and X 2 (t, s) are the same for both the Q and O 
operators (one can easily check that they satisfy the same 
equations), and the difference comes from other terms. In 
the cases where X\(t, s) and X 2 (t,s) are dominant, one 
would not expects sharp difference between the fermionic 
and bosonic baths. For example, when wi = w%, two 
operators Q and O are exactly the same. However, we 
found that the extra terms X${t, s) and X±(t, s) occurred 
in O may become important under certain conditions as 
shown in Fig. @] where the coefficients in the O and Q 
operators are plotted. Clearly, the fermionic and bosonic 
baths may result in very different dynamics. In the nu- 
merical simulations, we choose the Ornstein-Uhlenbeck 
noise K(t,s) = 2 e -(7+*n)|t-s| f or snTL p nc ity. However, 
our approach is applicable for arbitrary kinds of correla- 
tion function. 
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FIG. 4: Dynamic evolution of the coefficients in Q (for 
fermionic bath) and O (for bosonic bath) operators. The pa- 
rameters are u>i — 2, W2 = w = 1, 7 = 0.4, f2 = tt/4. 



VIII. CONCLUSION 

In summarizing, we have developed a novel technique 
called the fermionic quantum state diffusion approach 
which is a useful tool for studying quantum open sys- 
tems coupled to a fermionic bath. Using the Grassmann 
coherent state, the exact fermionic QSD equation and 
the corresponding master equation are derived for sev- 
eral physically interesting models. We have shown that 
the time-local QSD approach developed in this paper 
can efficiently solve open systems coupled to fermionic 
baths by employing the exact or approximate Q opera- 
tors. Moreover, our research also suggests that some spin 
bath problems can also be solved by using the effective 
fcrmion bath. Finally, it is of great interest to apply the 
fermionic QSD approach to more realistic models such 
as finite temperature fcrmion baths and large spin baths, 
and we leave these topics open for future discussion. 
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Appendix A: Spin chain as an effective fermion bath 

model 

In this section, we consider a spin-chain model where 
some spins are treated as the system of interest, the rest 



is treated as its environment. We show that the model 
can be transformed to a fermionic bath model. 

Consider a quantum system interacting with a XX spin 
chain. The Hamiltonian is 

H tot = H s + H h + H lnt , (Al) 

#b = ]>>+ar + i+^Vn, (A2) 

i 

H int = Da x +cr\l. (A3) 



After performing the Jordan- Wigner transformation, 

•,-, (A4) 



j'-i 
aj = exp(-J7ry^c^Cfc)i 



'i> i 



and the Fourier transformation 26] , 

N/2 

a j = TXf H exp(-ij4> p )a t 

ViV p=-N/2 

the original Hamiltonian Eq. (jAHIA3[) become 

Htot — H s + Hh + Hint , 

N/2 



(A5) 

(A6) 

(A7) 



H\y = 2_^ 2cos(f> p a,p dp, 

p=-N/2 
N/2 

Hi n t = —= 2_j [-£*exp(— i<j> P )a p + ex.p(i<t) p )a p L]. 

V N p= -N/2 

(A8) 
This effective Hamiltonian obtained from the transforma- 
tion takes the same form given by Eq. (JT][3|). Therefore, 
we may use the QSD approach to study the dynamics of 
the subsystem of the spin-chain model. 



Appendix B: Derivation of the non-Markovian QSD 
equation for a fermionic bath 

To start with, we list several useful commutation rela- 
tions between fermionic coherent state and operators: 

(m = -^U &|aj = &l£ = £&|. (bi) 

Using these relations, we can derive the QSD equation as 
d 



Of 



MC) 



= -i<£|#tot(*Mot0)) 

= -i{£\H 8 + J2(9te iUit 4L + H.C.)\Aot(t)) 






(B2) 



10 



where Q — —1^2^*^^^*. Then using the chain rule to 
introduce the functional derivative, 



m\Tptot(t)) = ^Mn 



(B3) 

(B4) 



Finally, we have 



§;MC) = [-i&. + kt - V J dsK(t, s ) A]^ (r); 

(B5) 
where K(t,s) = EN 2 e~ 8a,i(t_s) - This is just the final 

i 

QSD equation. 

Appendix C: Equation for Q operator 

First, consider the following two commutation rela- 
tions: 

^:mc)} = -^-^mo (ci) 

and 

§zWt(n] = QhOQMO + (^Q)MO, (C2) 

for fixed order of Q operator. One can prove them easily. 
With Eq. (jCll) and Eq. (jC2[l . we can apply the con- 
sistency condition to ipt(£,*)- 

£!*«*> = I4* ( °- (C3) 



The left-hand side is 



LHS = [-iH s Q - LC t 



-L^—Q)-DQ{-C)Q]Mn- (C4) 

On the other hand, the right-hand side becomes 
d 



RHS = -(Q)^ t (C) + {-iQH s + QLQ ~ QL^Q}MO- 

(C5) 
Equate LHS and RHS and eliminate ipt(£,*), we have 



^- t Q = -i[H s ,Q]-{L$,Q} 



Si 



VQ{-OQ + Ql ] Q-l^q. (ce) 



Appendix D: Proof of Novikov-type theorem for the 
Grassmann noise 

In this section, we will provide a proof of Novikov-type 
theorem for a Grassmann Gaussian noise, which plays a 
crucial role in deriving the exact or approximate master 
equations from the corresponding stochastic Schrddinger 
equations. 

Theorem: Suppose that £t,£* are Grassmann-type 
Gaussian processes and the Pt is the stochastic density 
operator, then we have the following two identities: 



(PA} s = -(Q(C)Pt)s, 
($P t ) g = (P t QH-0)s- 



Proof: 



(Dl) 
(D2) 



(P&h 



= JUd&dtie-Xi&i |V,(D> (^(-01 (*£ 



9j< 



%) 



Y.^-^JU^dtme)) w-o\ ^-^ aii 



di 



d(-a) 



MO) W-OIK 



j J * 

-* E a^/lL*?^ [«" E > ^ (/ *|§ 4 } A] 
-jdaY, \9i I 2 ^ wfUdGdti [e- a *6 |A] 



•Eitffc 



T[.d&dZi 



•EiCS« 



dsK{t,s)Q(t,s,e)Pt. 



■{QPth 



(D3) 
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Similarly, we can prove 

<£*p t ) s = (PtQH-Z))* 



(D4) 



This concludes our proof of the Novikov-type theorem for 
the Grassmann Gaussian noise. 



Appendix E: Exact Q operator for the two-qubit 
model 



where the basis operators are given by 



Q\ = a_, Q 2 = a_, Q 3 = 

A A-B -A A o^A~B 



^A±.B 
<J Z <7_, 



(E2) 



For the coupled two-qubit model, the exact Q takes 
the following form: 

Q(t,s,C) = A{t,s)Qi + f2(t,s)Q 2 + h(t,s)Q 3 



and fj (j — 1,2,3,4,5) are some time-dependent coeffi- 
cients. Substituting Eq. (IE1|) into Eq. (fTTj) . we obtain a 



+ f 4 (t,s)Q 4 + i 



I ds' /k(i s s')£*iQ5 (El) se t of partial differential equations governing the coeffi- 
Jo cients of the Q operator, 



q^/i(*> s ) = + 2iuj Afi + n A F 1 f 1 + n B F 3 f l - iJ xy f 3 - n B F 1 f 3 + k b F 4 / 3 + 2iJ z f 4 + n A F 4 f 4 + n B F 3 f 4 - in B F 5 , 

d 

o^/aC*) s ) = + 2iuj b/ 2 + K A F 4 f 2 + K B F 2 ,f 2 + liJzh + k a F 4 J 3 + n B F 3 f 3 - iJ xy f 4 ~ n A F 2 f 4 + n A F 3 f 4 - ik a F 5 , 

d 

ol/sv*) s ) = - % Jxyfi - K A F 2 f! + K A F 3 fl + 2lJ z f 2 + K A F 4 f 2 + K B F 3 f 2 + 2lLU B f 3 + K A F 4 f 3 + K B F 2 f 3 - IK A F 5 , 

d 

fllfa(t> s ) = +2lJzfl + K A F 4 f 1 + K B F 3 f 4 - lJ xy f 2 - K B Flf 2 + K B F 4 f 2 + 2lLU A f 4 + K A F 4 f 4 + K B F 3 f 4 - IK B F 5 , 

d 
— f 5 (t, s, S ) = +K A F 5 fi + K B F 5 f 2 - K B F 5 f 3 - K A F 5 f 4 + 2iuj A f 5 + 2iuj B f 5 + K A F 1 f 5 + n A F 4 f b + K B F 2 f 5 + K B F 3 f 5 , 

(E3) 



where F (t) = £ dsK(t,s)fj(t,s) {j = 1,2,3,4) and 

Fs(t, s') = LdsK(t,s)f5(t,s,s'), with the initial con- 
ditions: 

fl(t,8 = t) = K.4, 
f 2 (t,S = t)= K B , 
h(t,8 = t) = 0, 
U(t,8 = t) = Q, 

f 5 (t,s = t,s') = 0, 

f 5 (t,S,s' = t)=i[K A f2(t,s)+KBh(t,s)]. (E4) 



Appendix F: Differential equations for coefficients of 
bosonic O in example 4 

The coefficients in Eq. (1851) satisfy the following dif- 
ferential equations 

d 

-X 4 (t, s) — iuJ a Xl +X1X1 + x 2 Xi, (Fl) 



dt 
di 



d_ 
di 



d_ 
di 



xa(t, s) = iujbx 3 — x 4 X 2 + x 3 X 2 + x 2 X 3 
+ x 3 X 3 - x 3 X 4 - x 2 X 4 - iX 5 , 



x 4 (t, s) = iui a x 4 + x 4 X\ + x\X 4 — x\X 3 
— x 3 X\ + x 4 X 3 — x 4 X 4 — 1X5, 



0_ 

dt 



x 5 (t, s, s') = iuj a X5 + iu)bX5 + X5X1 + x 5 X 2 



(F3) 



(F4) 



x 2 {t,s) = iuj b x 2 + xiX 2 + x 2 X 2 , 



(F2) 



+ x 5 X 3 - x 5 X 4 + x 4 X 5 + x 2 X 5 , (F5) 

with the initial conditions 

x 1 {t,s = t) = 1, (F6) 

x 2 (t,s = t) = l, (F7) 

x 3 (t,s = t)=0, (F8) 

x 4 (t,s = t) = 0, (F9) 

x 5 (t,s = t,s') = 0, (F10) 

ix5(t, s, s' = t) = 2(x 2 — x\) + x 3 + x 4 . (Fll) 



and 
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Xj(t) = K(t,s)xj(t,s)ds (j = 1 to 4) (F12) 
Jo 

X 5 (t,s')= f K(t,s)x 5 (t,s,s')ds (F13) 

Jo 
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